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Abstract 

The resolution of Dyson-Schwinger equations leads to the freezing of the QCD run- 
ning coupling (effective charge) in the infrared, which is best understood as a dynamical 
generation of a gluon mass function, giving rise to a momentum dependence which is 
free from infrared divergences. We calculate the interquark static potential for heavy 
mesons by assuming that it is given by a massive One Gluon Exchange interaction 
and compare with phenomenologyical fits inspired by lattice QCD. We apply these 
potential forms to the description of quarkonia and conclude that, even though some 
aspects of the confinement mechanism are absent in the Dyson-Schwinger formalism, 
the spectrum can be reasonably reproduced. We discuss possible explanations for this 
outcome. 
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1 Introduction 



The discovery of asymptotically free constituents of hadronic matter in deep inelastic 
scattering experiments gave birth to Quantum Chromodynamics (QCD) pQ as the 
accepted theory of strong interactions. However, the most precise experimental data, 
the hadron spectrum and hadron properties, have been obtained in the low-energy 
region. While QCD explains asymptotic freedom [2j|3], it has not been shown to 
describe confinement of quarks and gluons, nor the realization of chiral symmetry. 
Therefore, predicting low-energy properties of strongly interacting matter represents a 
theoretical challenge. 

The development of nonperturbative techniques are essential for undertsanding 
QCD in this regime. Lattice gauge theory [HE] constitutes a nonperturbative reg- 
ularisation scheme which allows numerical solutions of the theory describing properties 
of interacting QCD matter The accuracy of lattice results has been tremendously 
improved during the past decade with the availability of more powerful computers 018], 
and lattice results are considered in many instances the data which other nonpertur- 
bative schemes should reproduce. 

The approximate resolution of the Dyson- Schwinger (DS) equations is another non- 
perturbative approach which has progressed considerably in the last ten years in great 
part due to the interplay between their findings and lattice results. It is a more ana- 
lytical approach and has led to a very appealing physical picture establishing that the 
QCD running coupling (effective charge) freezes in the deep infrared. This property can 
be best understood from the point of view of a dynamical gluon mass generation [9|[T0] . 

The aim of this presentation is to investigate the consequences associated to the 
static interaction one can derive from this picture. For this purpose we calculate numer- 
ically the one gluon exchange (OGE) static potential deriving from the DS equations 
and we compare it to phenomenological potentials whose shape has been inspired by 
lattice computations. The application of these potentials to the description of quarko- 
nia is discussed. 

The paper is written as follows, in the next section we describe some phenomeno- 
logical potentials justified or motivated by lattice calculations. In Sec. 3 we obtain a 
potential from the resolution of DS equations and describe some of its characteristics. 
In Sec. 4 we make a comparative analysis of the different potentials in their application 
to quarkonia. Finally, in Sec. 5, we summarize our main findings and discuss a possible 
interpretation of our results. 

2 Heavy quark dynamics from Lattice 

Soon after the discovery of the J/\l/ mesons in e + e~ annihilation, the possibility of a 
nonrelativistic treatment of such states, in analogy to the positronium of electrody- 
namics, was suggested [11] . Quarkonia, i.e. mesonic states that contain two heavy 
constituent quarks, either charm or bottom, owe their name to this analogy. 
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Figure 1: We show the three potentials described in the text : Cornell (solid), Aachen 
(dashed) and SESAM (dotted). For the Cornell and Aachen potentials the parameters 
have been chosen phenomenologically to provide a reasonable description of the spec- 
trum and their values are given in Sec. 4. The parameters of the SESAM potential are 
chosen to reproduce the shape of the Cornell potential at short distances. 



For sufficiently heavy quarks, it can be shown in lattice QCD, that the bound state 
problem becomes essentially nonrelativistic and the dynamics is controlled approxi- 
mately by a Schrodinger equation with a static potential [12j[T3]. Lattice calculations, 
without dynamical quarks (quenched), give rise to a static potential of the form [T2"|[T3] 

V(r) = —a/r + br. (1) 

containing the perturbative expected Coulomb term plus an additional linear term. 

The Cornell group, prior to the lattice QCD derivation, using a and b as parameters, 
successfully applied this potential to the phenomenological description of the low lying 
quarkonia states [141417] and therefore this potential function is known as the Cornell 
potential. 

Another static interquark potential shape has been guessed by performing un- 
quenched lattice QCD calculations. Using Kogut-Susskind fermions and adopting the 
lattice spacing from the p mass measurements a potential with the following structure, 

V( r ) = {-a/r + 0r) (— ^— ) , (2) 

has been shown to describe these lattice results [18J. This parameterization shows 
screening, associated to quarks loops, which has been cast in the form of the additional 
factor multiplying the Cornell potential, leading at large values of r to a constant 
P/j. Henceforth we shall call this shape the Aachen potential and we shall choose its 
parameters to provide a reasonable fit to the spectrum. 
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More recent unquenched lattice calculations have been performed by the SESAM 
collaboration [19]. The mixing and the transition (string breaking) between qq and 
MM (M representing a meson) components have been studied. A correlation matrix 
for this two-level problem has been considered, the matrix elements incorporating light 
quark propagators into the standard Wilson loop. Compelling evidence, both for ex- 
plicit mixing and string breaking, has been demonstrated so that for r > r c > 1 fm (r c 
representing the string breaking distance) the MM configuration becomes energetically 
favorable. The static qq potential exhibits screening and saturates at twice the mass 
of the M meson. An approximate parametrization of it for r < r c has been proposed, 

V(r)=m M + m JI + g(r)G(r) + C (3) 

where 



g(r) = arctan [d(r — r c )] 




with C, r c , d, e and a obtained from lattice results. This parameterization has not the 
correct long distance limit (rriM +^m)' anc ^ differs little from the Cornell potential for 
r < 0.75 fm. 

We show the three potentials in Fig. [TJ for comparison^. We use in what follows 
the first two, which contain the two aspects we want to emphasize in here, linear 
confinement and screening. 



3 Heavy quark dynamics from Dyson-Schwinger Equa- 
tions 

Infrared finite solutions for the gluon propagator of quenched QCD are obtained from 
the gauge-invariant nonlinear Dyson-Schwinger equations formulated in the Feynman 
gauge of the background field method. These solutions may be fitted using a massive 
propagator J9j[10]. Even though the gluon is massless at the level of the fundamental 
QCD Lagrangian, and remains massless to all order in perturbation theory, the nonper- 
turbative QCD dynamics generates an effective, momentum-dependent mass, without 
affecting the local SU(3) C invariance, which remains intact. 

The gluon mass generation is a purely nonperturbative effect associated with the 
existence of infrared finite solutions for the gluon propagator, A(g 2 ), i.e. solutions with 
A _1 (0) > 0. Such solutions may be fitted by a "massive" euclidean propagator of the 
form PEI] 

1 The spin-dependent corrections to the potentials above can also be derived from lattice QCD [20 . 
We do not consider them here. 
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A(g 2 ) 



1 



(4) 



■ 2 + m 2 (q 2 ) ' 



where m 2 (q 2 ) depends nontrivially on the momentum transfer q 2 . One physically moti- 
vated possibility, which we shall use in here, is the so called logarithmic mass running, 
which is defined by 



where m ,p and 5 are parameters whose values are chosen to fit the lattice propagator 
and A is the QCD scale. Note that in the limit q 2 — > one obtains m 2 (0) = mf), giving 
meaning to m®. 

Because of the presence of this dynamical gluon mass the strong effective charge 
extracted from these solutions freezes at a finite value, giving rise to an infrared fixed 
point for QCD (9HTU] - This nonperturbative generalization of ct{q 2 ), the QCD running 
coupling, comes in the form 



where a = f- and we take (3 = 11 — 2n//3 where nf is the number of flavors. Note that 
its zero gluon mass limit leads to the LO perturbative coupling constant momentum 
dependence. The m(q 2 ) in the argument of the logarithm tames the Landau pole, and 
a(q 2 ) freezes at a finite value in the IR, namely a _1 (0) = /3ohi(pm 2 (0)/A 2 ) . 

Let us construct a simple potential model where the main source of dynamics is the 
One Gluon Exchange potential (see Fig. [2]) with the propagator and coupling defined 
by Eqs. EJE]andEJ 



The potential between static charges is related to the Fourier transform of the 
time-time component of the full gluon propagator as, 







Figure 2: One gluon exchange contribution to the potential. 




(7) 



4 



where Cp is the Casimir eigenvalue of the fundamental representation of SU(3) [Cf — 
4/3], the bold terms, q and r, are 3- vectors and A o(q) is the zero-zero component of 
the gluon propagator in the momentum configuration. 

Using the above equations the expression a(q 2 )A o(q) in the integrand of EqfT] 
becomes 



d(q 2 



q 2 + m 2 (q 2 ) 



with the dynamical mass m 2 (q 2 ) and the nonperturbative coupling constant a(q 2 ) 
defined by Eqs. [5] and [6] . This particular combination of propagator and coupling 
constant is renormalization group invariant. The structure of <i(q 2 ) is shown in Fig. |3j 




q/A 



Figure 3: Numerical solution for the RG invariant product d for m = 360 MeV, p — 1, 
A = 300 MeV and rif = 4. We plot the adimensional function dA 2 as a function of 
q/A. 



In spherical coordinates EqJT] becomes, 
V(t)= 1 ' 



d(j) / rf|q|q 2 a(q 2 ) A(q 2 ) / d(cos6)e 

7T JO JO J-1 



i|q| |r| cos ( 



By performing the angular integral we obtain, 



V(r) 



8C f 



r Jo 



d\q\ |q| a(q 2 ) A(q 2 ) sin(|q| |r|). 



(9) 



(10) 



The complete QCD potential involves multigluon exchanges, but we assume that 
the leading nonperturbative term will be sufficient for the description of the spectrum. 

In Fig. H]we show the potential (up to a constant) derived from Eq. [TDl i.e. from 
the DS equations with the definitions in Eqs. [5] and [61 We have chosen the following 
range of parameters: m ~ 360 - 480 MeV, p = 1 - 4, 5 = 1/11 0[22l[23]. The value 
of A has been taken to be 300 MeV. In order to adjust the behavior at the origin to 
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Figure 4: We show the potential obtained from the DS equations. The shaded range 
corresponds to m ~ 360 — 480MeV and p ~ 1 — 4. The figure on the left corresponds 
to 4 flavors( (3 = 25/3). The figure on the right corresponds to 5 flavors (/3 = 23/3). 
For comparison we plot the Cornell (dashed) and Aachen (dotted) potentials. 



the data we have used /3q corresponding to 4 and 5 flavors. To do this appropriately 
one should introduce the running of the quark masses, which are at present not well 
known. However, since asymptotically the masses run to zero, our way of proceeding 
achieves the correct strength of the potential at low r. The potential describes well 
the low radial behavior, by construction, and flattens at large r going asymptotically 
to zero. 

We should remove from the above calculation the additive infinite self-energy contri- 
bution associated with the static sources [12]. In lattice QCD this is done normalizing 
the potential such that V(tq) = where r$ is the Sommer scale [24] . We proceed 
in the same way but take the subtraction point at the zero of the phenomenological 
potentials, which happens to be at r$ ~ 0.35/m (see next section). The result of this 
procedure is shown in Fig. [5j In this way we increase the value of the potential without 
changing its shape. We call the result of this construction the Dyson- Schwinger (DS) 
potential . We obtain in this way a potential which resembles the Aachen and not the 
Cornell potential. 

It should be noted that the shape of the DS potential does not vary considerably 
when we change the parameters within the expected theoretical range. There is no way 
to reproduce the large r behavior of the Cornell potential, instead the DS potential 
flattens and becomes asymptotically constant, similarly to the screened potentials. 
However, it should not! The approximations used to find the solution to the DS 
equations do not contain quark loops and therefore they incorporate no mechanism 
for screening, i.e. a mechanism derivable from the breaking of the string [18,25j. The 
truncation in the DS set of equations and the absence of multigluon exchanges might 
be responsible for the missing linear rise at large r. 
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Figure 5: We show the DS potential after the Sommer subtraction, for n/ = 4 and the 
same range of parameters as in Fig. 0] For comparison we plot the Cornell (dashed) 
and Aachen (dotted) potentials. 

4 Quarkonia description by static potentials 

The Cornell potential was introduced for phenomenological reasons predating the com- 
plete lattice derivation. From the spin-averaged quarkonia spectra it was evident that 
the underlying potential could not be purely Coulomb type [H] . Therefore, the po- 
tential was implemented by a sum of the perturbative expectation plus an additional 
linear term, recall Eq. [TJ Aiming at a universal treatment for charmonium (cc) and 
bottomonium (bo) the Cornell potential has been used with the same values of the 
parameters, a and b, in both cases, for the quark masses, m c and ra&, respectively. A 
typical range of values providing a reasonable fit to the masses of the low lying states 
(up to 1.0 GeV excitation energy) is a ~ 0.51-0.52 and \fb ~ 412-427 MeV (I5HI7]. 
An example of such spectral fit for charmonium is provided in Table 1 where we have 
chosen \[b = 427 MeV and a = 0.52 with m c = 1350 MeV and compared the results 
with masses of experimental resonances having well established J PC quantum numbers, 
most of them with J PC = 1 coming from ISR (Initial State Radiation) processes. 
The Cornell model provides a good fit to the lower states (at most 30 MeV differ- 
ence with data) but it cannot accommodate all the known higher energy resonances 
but only some of them (we use charmonium instead of bottomonium to clarify the ef- 
fect). For instance ?/>(4040), ?/>(4160) and ■0(4415) maybe assigned to the 3s, 2d and 4s 
states respectively. Then other two resonances, cataloged in the Particle Data Group 
Review (3Tj as A (4260) and A (4360), cannot be fitted. 

The screened potentials introduced initially to fit lattice data have been used phe- 
nomenologically [28j[29]. However, the applicability of screened potentials to the spec- 
tral description has been a matter of debate [30] since their use is not justified above 
the meson-meson string breaking threshold. We assume that in the case of quarkonia 
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3502 
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3525.3 ±0.2 
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2p 


3983 


3894 
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Table 1: Calculated charmonium masses, M Co meih M Aachen and M DS from the Cornell, 
Aachen and DS potentials. For Cornell and Aachen a = a = 0.52 andVb = V^ = 427 
MeV. The remaining parameter in the Aachen potential has been chosen to be 7 = 0.38 
fm _1 . For DS p = 1 and mo = 345.7 MeV. The charm masses are m c = 1350 MeV for 
Cornell and m c = 1400 MeV for Aachen and DS. Masses for experimental candidates, 
Mp£,Gi have been taken from [31]. For p waves we quote the centroid of the npo, npi 
and np 2 states. 



there is an effective string breaking threshold sufficiently high in energy, to allow for 
a description of the known spectrum in terms of screened potentials. For the sake of 
comparison the results from the Aachen potential for m c = 1400 MeV and the same 
values for the corresponding Cornell parameters a = a, (3 = b are presented in Table 
1. The value of 7 = 0.38 fm _1 is fixed to get a reasonable spectral fit. It is important 
to emphasize that the values of the parameters of the Aachen potential extracted from 
lattice data [18] can only give a reasonable description of the masses of the low lying 
quarkonia. 

One should keep in mind that these potentials do not contain spin-dependent terms 
which makes them reliable only when these terms do not play a major role. We use 
here the conventional approximations, e.g. we consider that these potentials should fit 
better the spin triplet states and take the centroids of p states as data for comparison 
with our results. 

The main difference between the Cornell and Aachen potentials refers to the de- 
scription of the higher excited states. The Aachen potential may allow a one to one 
assignment of the calculated states to the data. 

The similarity of the DS potential (for certain parameter sets) to the Aachen po- 
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tential as shown in Fig. [5] motivates the exercise of fitting the spectrum with the DS 
potential. As can bee seen in Table 1 this can be achieved for a set of parameters which 
is close to the expected theoretical range. A low value of the mass, mo = 345.7 MeV, 
is necessary since only with a value close to A one gets sufficient strength to achieve, 
after Sommer subtraction, an asymptotic behavior close to the Aachen potential used 
in our spectral fit. 

Let us discuss the limitations of the parameters used in the fit. The 0q expression 
is fixed by QCD. /3 is the leading order coefficient of the beta function of the theory, 
which is scheme independent; the small difference in value implied by the choice of n/ 
has almost no effect on the results. A is a scale fixed by dimensional transmutation. We 
take 300 MeV and keep it fixed. The values obtained in experimental fits range from 
250 — 300 MeV but we have not even used that liberty. The 5 value and the variations 
in the values of p and m are within the limits provided by the DS calculations, in 
particular A < m < 2A [10,21]. 

The fit of the spectrum is strongly dependent on the detailed shape of the potential 
as determined by the mass function and coupling constant, Eqs. [5] and [6] respectively. 
Recall that the mass function is a fitting function to the lattice results for the prop- 
agator [21]. The quality of the spectral fit seems to point out, that the DS potential 
may qualitatively be a good description of the dynamics. 

5 Discussion and Conclusions 

In order to interpret the above results we have to resort again to lattice calculations. In 
ref. [2B] the authors investigated the origin of the long range linearly rising potential. 
Using a mechanism described in ref. [27] they were able to subtract the contribution 
of the center vortices leading to a flat potential as shown in Fig. [6] on the left. On the 
right of Fig. [6] we show the DS potential after Sommer subtraction and the linear rise 
of the Cornell potential. It is quite apparent that the figure resembles the one on the 
left. Therefore we may tentatively conclude that the DS potential contains the physics 
associated with the approximate gluon interaction but does not contain the physics of 
the confinement mechanism. We do not advocate any mechanism for confinement in 
QCD, we only adhere to the fact that the confinement mechanism, whatever it be, is 
the mechanism behind the rising potential in the quenched approximation. 

A surprising result of our calculation is the actual similarity of the DS potential 
to the Aachen potential as shown in Fig. [5J The latter arises due to the breaking of 
the string and is represented by a screened potential [IB] . The similarity is astonishing 
more so since we have used conventional values for p = 1,8 — 1/11, A = 300 MeV and 
only varied mo, which becomes for the best fit mo = 345.7 MeV within the expected 
range A < m < 2A [T01I2T] . 

Not aiming at such precision if we fit the upper extreme of the shaded region, where 
the DS potential is defined by mo = 360 MeV and p = 1, to the Aachen potential shape 
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Figure 6: The figure on the left corresponds to the calculation of ref. [26]. In the 
figure on the right the dotted line corresponds to the DS potential for n/ = 4 and 
mo = 360 MeV and p — 1. The dashed line corresponds to the linear piece of the 
Cornell potential. The full Cornell potential is drawn (solid line) for comparison. 



function we get a = 0.40, a//3 = 412 MeV and 7 = 0.53 fm _1 . The other extreme set of 
parameters for the DS potential, the lower edge of the shaded region, m = 480 MeV 
and p = 4, leads to a = 0.37, = 361 MeV and 7 = 1.81 fm _1 . It is clear then, 
that the main difference between the parametrizations arises due to the value of the 
different screening ranges, which in the Aachen potential is controlled by 7. 

Is this similarity accidental or does string breaking imply a dilution of the confine- 
ment mechanism associated with multigluon exchanges and/or higher order truncation 
schemes ? Does the Sommer subtraction introduce the energy scale of string breaking 
? More research needs to be done in order to understand the confinement and string 
breaking mechanisms though our investigation hints a possible scenario. 

In conclusion, we have calculated the OGE potential associated to the approximate 
resolution of the Dyson- Schwinger equations for the gluon propagator. The low r be- 
havior is determined perturbatively. The large r behavior is certainly nonperturbative. 
The Sommer procedure, to avoid self-energy effects of the static charges, leads to a 
potential which is not negative at large r. The approximate resolution of the DS equa- 
tions is not able to reproduce, at the level of the OGE, the Cornell potential. The DS 
equations together with the Sommer subtraction is close to the Aachen potential that 
contains a string breaking mechanism. The fact that we have to push the parameters 
to the limit of the allowed region to reproduce the spectrum maybe due to the fact 
that we have not treated the quark mass terms appropriately. 

The potential derived from the DS equations might contain most of the dynamics 
associated with the interquark interaction apart from a non trivial constant which 
might be related, in the unquenched calculation, to the confinement mechanism. 
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